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ABSTRACT 


Dissipation  rates  of  turbulence  kinetic  energy  e  and  enstrophy^  are  reported  in  a  high 
Reynolds  number  turbulent  wake.  Previous  turbulent  wake  observations  have  been  made 
in  laboratory  experiments  with  relatively  low  Reynolds  number  flows,  O(103).  Our  results 
represent  a  set  of  rare  field  observations  of  vorticity  and  turbulence  in  a  turbulent  wake 
with  a  high  Reynolds  number,  0(1O7).  The  turbulent  wake  was  formed  by  an  unsteady 
strong  tidal  current  interacting  with  a  bridge  pier.  Measurements  were  taken  in  the 
intermediate  wake.  Both  e  and  £  show  a  similar  downstream  decay  rate  faster  than  that 
predicted  by  the  self-preservation  similarity  in  the  far  wake.  The  theoretical  relation 
e  -  v£2  for  high  Reynolds  number  flow  is  confirmed  by  field  observations.  The 
magnitudes  of  the  vertical  and  horizontal  components  of  enstrophy  are  not  significantly 
different.  Probability  density  functions  of  normalized  vorticity  components  and  of  the 
dissipation  rate  of  turbulence  kinetic  energy  are  nearly  identical  to  each  other.  They  are 
also  nearly  identical  to  the  empirical  function  found  in  laboratory  experiments.  The 
turbulence  internal  intermittency  is  ~0.2,  estimated  from  autocorrelation  coefficients  of 
enstrophy;  this  value  is  close  to  that  reported  previously  in  turbulent  wakes  and  jets. 
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1.  Introduction 


Turbulent  wakes  have  been  studied  extensively  in  laboratory  experiments, 
numerical  models,  and  theoretical  studies.  Field  observations  of  turbulent  wakes  are 
difficult  to  obtain  because  the  natural  environment  is  not  controllable,  the  spatial  scale  is 
too  large  to  cover  synoptically,  and  the  turbulent  wake  is  often  unsteady.  Natural  turbulent 
wakes  often  have  a  Reynolds  number  Rej  =ux  d  u"1,  where  m*  is  the  ffee-stream  velocity,  d 
the  scale  of  the  obstacle,  and  u  the  kinematic  molecular  viscosity,  many  decades  greater 
than  those  laboratory  and  numerical  experiments  can  achieve.  Whether  results  obtained  in 
the  relatively  low  Reynolds  number  flows  are  applicable  to  high  Reynolds  number  natural 
flows  remains  a  question. 

In  turbulent  wakes  vorticity  plays  the  most  important  dynamic  role.  The  vorticity  is 
generated  via  the  flow  field  interacting  with  the  obstacle  and  cascades  to  larger  scales. 
Unfortunately,  vorticity  is  a  difficult  variable  to  obtain  because  it  requires  accurate 
measurements  of  velocity  gradients.  In  laboratory  experiments  Antonia  et  al.  [1]  designed 
a  probe  to  measure  vorticity  in  the  wake  of  a  circular  cylinder  with  a  Reynolds  number 
0(  1 03).  Sanford  et  al.  [2]  presented  the  first  set  of  oceanic  vorticity  measurements  in  a  tidal 
channel  with  a  Reynolds  number  O(107),  using  a  unique  electro-magnetic  vorticity  meter 
(EMVM)  that  also  measures  velocity  and  dissipation  rate  of  turbulence  kinetic  energy  e. 
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Here,  we  present  observations  of  vorticity^  and  turbulence  kinetic  energy 
dissipation  rate  e  in  a  turbulent  wake  behind  the  pier  of  the  Tacoma  Narrows  Bridge  in 

n 

Washington  State,  with  Red  >2x10  .  The  turbulent  wake  was  established  by  a  strong, 
unsteady,  decelerating  tidal  current  in  the  presence  of  a  weak  cross-wake  shear. 

The  experiment  and  measurements  are  described  in  section  2.  General  statistics  of 
vorticity  and  turbulence  kinetic  energy  dissipation  rates  are  presented  in  section  3.  The 
downstream  variation  of  £  and  self-preservation  properties  are  discussed  in  section  4.  The 
downstream  variation  of  ens trophy^  and  turbulence  intermittency  determined  from  their 
autocorrelation  coefficients  are  described  in  section  5,  and  a  summary  follows  in  section  6. 

2.  Experimental  Conditions 

A.  Instruments  and  Measurements 

Measurements  of  velocity,  vorticity,  and  turbulence  kinetic  energy  dissipation  rates 
were  taken  behind  a  pier  of  the  Tacoma  Narrows  Bridge  (Fig.  la)  during  a  brief  study  in 
southern  Puget  Sound.  Within  3  h,  three  sets  of  zigzag-transect  observations  were  taken 
downstream  from  the  pier  during  the  decelerating  phase  of  the  northward  ebb  tidal  current. 
Tidal  currents  in  Tacoma  Narrows  often  exceed  2  ms'1.  An  aerial  photo  of  the  old  Tacoma 
Narrows  Bridge,  taken  >50  years  ago  soon  after  its  failure,  reveals  vortices  generated 
behind  the  pier  (Fig.  lb).  Similar  vortices  and  turbulence  formed  in  the  wake  of  the  new 
Tacoma  Narrows  Bridge  pier  are  the  focus  of  this  analysis. 
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Our  primary  instrument  is  the  EMVM  tow  body  (Fig.  lc).  Sensors  on  the  EMVM 
include  two  EM  velocity  and  vorticity  meters  (VMs),  a  microstructure  shear  probe,  and  a 
CTD  sensor.  Details  of  sensors  and  measurements  are  described  in  Table  1 .  The  principle 
of  measuring  vorticity  and  velocity  using  the  VM,  and  sensor  response  functions  are 
described  by  Sanford  et  al.  [2]. 

One  VM  is  mounted  on  the  bottom  of  the  EMVM  (Fig.  lc;  labeled  EM  velocity  and 
vorticity  sensor).  It  measures  the  streamwise  (i.e.,  along  the  axis  of  the  EMVM)  and  the 
spanwise  (transverse  to  the  axis  of  the  EMVM)  components  of  velocity  and  the  vertical 
component  of  vorticity.  Another  VM  is  mounted  on  the  side  of  the  EMVM.  It  measures  the 
streamwise  and  the  vertical  components  of  velocity,  and  the  spanwise  component  of 
vorticity.  Measurements  of  s  are  obtained  from  the  microstructure  shear  probe  in  front  of 
the  EMVM.  Ancillary  sensors  measure  the  pitch,  compass,  pitch  rate,  and  yaw  rate.  These 
ancillary  measurements  are  needed  to  remove  the  instrument  motion  from  velocity  and 
vorticity  measurements. 

Because  of  the  finite  size  of  the  vorticity  sensor,  variations  of  vorticity  and  velocity 
at  scales  smaller  than  the  sensor  separation,  0.09  m,  are  attenuated  in  VM  measurements 
[2],  The  turbulence  velocity  has  red  spectral  characteristics,  i.e.,  large  eddies  contain  most 
of  the  energy.  Turbulence  velocity  measurements  are  less  affected  by  the  sensor 
attenuation  at  small  scales.  The  turbulence  vorticity  has  a  blue  spectrum.  Most  of  the 
enstrophy  is  at  scales  close  to  the  Kolmogorov  scale,  rj  =  v3/4e~,/4.  Because  the 
Kolmogorov  scale  decreases  with  increasing  turbulence  intensity,  the  amount  of  sensor 
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attenuation  on  our  vorticity  measurements  increases  with  turbulence  intensity.  Sanford  et 
al.  [2]  derive  response  functions  of  the  sensor  attenuation  for  velocity  and  vorticity.  They 
demonstrate  that  the  sensor  attenuation  effect  can  be  corrected.  In  the  following  analysis, 
response  functions  of  the  vorticity  and  velocity  are  applied  to  our  measurements  to  correct 
for  the  missing  variances  at  small  scales. 

The  EMVM  was  towed  at  speeds  of  1-2  m  s-1  through  the  water  at  a  depth  of  ~5.5 
m  in  water  of  1 5-40  m  depth  (Fig.  la).  Measurements  were  taken  from  upstream  of  the  pier 
to  2-3  km  downstream  of  the  pier.  The  pier  is  located  in  water  of  ~40  m  depth  and  has  a 
nearly  square  shape  of  ~30  m  on  each  side  (d  =  30  m).  Most  previous  wake  studies  have 
focused  on  turbulent  wakes  behind  a  cylinder. 

In  this  analysis  we  rotate  the  coordinate  system  into  the  along-wake  and 
across-wake  components  (Fig.  la).  The  along-wake  coordinate  is  defined  as  61  degrees 
counterclockwise  from  east,  and  the  across-wake  coordinate  is  90  degrees 
counterclockwise  from  the  along-wake  coordinate.  Similarly,  measured  horizontal  velocity 
is  rotated  to  the  along-wake  and  across-wake  components. 

B.  Free-stream  Velocity 

Our  measurements  were  taken  for  ~3  h  beginning  from  the  peak  of  an  ebb  tidal 
current.  To  obtain  the  estimates  of  the  free-stream  velocity  ux ,  we  compute  averages  of  the 
streamwise  velocity  taken  during  the  zigzag  transects  outside  of  the  wake  (Fig.  2).  The 
free-stream  velocity  decreased  from  2.5  m  s~*  to  0.5  m  s  1  within  3  h  and  exhibited  a  slight 
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transverse  gradient,  likely  due  to  the  lateral  boundary  friction.  We  fit  the  observed 
free-stream  velocity  to  a  model  consisting  of  a  constant  transverse  gradient,  a  mean  current, 
and  a  semidiurnal  tidal  current  described  as 


*C(>V)=  [1  -  ydyuj[u  +  u,  cos(<y/)],  (1) 

where  dyum  =  6.3  x  10-4  s~]  is  the  fitted  transverse  gradient  of  the  free-stream  velocity,  u 
=  1.09  m  s'1  is  the  fitted  mean  current,  u  =  1.32  m  s_1  is  the  fitted  tidal  current  amplitude, 
and  to  =  2rc/12.4  h  is  the  lunar  semidiurnal  frequency.  The  model-fit  free-stream  velocity 
<  is  used  to  estimate  the  turbulence  kinetic  energy  dissipation  rate  predicted  by  the 

turbulent  wake  self-preservation  similarity  and  to  compare  with  our  direct  observations  of 
turbulence  kinetic  energy  dissipation  rate.  During  the  three  lateral  transects,  the  ship  was 
set  downstream  by  the  tidal  flow  at  average  speeds  of  1.95,  1.87,  and  1.01  m  s-1, 
respectively.  The  average  free-stream  velocity  during  the  three  transects  was  2.2,  1.9,  and 
1 .0  m  s'1  through  the  water  (similar  to  the  ship’s  drift  velocity).  The  average  free-stream 
velocity  is  used  to  scale  observed  dissipation  rates  of  turbulence  kinetic  energy  and 
enstrophy  in  sections  4  and  5. 

3.  Statistics  of  Enstrophy  and  Turbulence  Kinetic  Energy  Dissipation  Rate 

In  high  Reynolds  number  flow,  the  dissipation  rate  of  turbulence  kinetic  energy  is 
related  to  the  turbulent  enstrophy  as3 
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where  C,1  is  the  total  enstrophy,  i.e.,  the  sum  of  the  squares  of  the  three  components  of 
vorticity.  This  relation  has  been  justified  in  low  Reynolds  number,  Re  =  1530,  flows  in  lab 
experiments  [4],  Our  measurements  were  taken  in  the  field  flow  with  Red  >2  x  10  ,  given 
Woo  >  1  m  s~\  d=  30  m,  and  y=1.4x  10~6  m2  s-1.  Therefore,  we  expect  (2)  to  be  satisfied. 

Only  one  of  the  two  components  of  horizontal  vorticity  was  measured.  We  assume 
that  the  vorticity  field  is  horizontally  homogeneous  and  estimate  the  total  enstrophy 
as£2  =  2 +  Cl ,  where  £%  is  the  spanwise  enstrophy  and  C,\  is  the  vertical  enstrophy. 

Large-scale  coherent  vortices  in  the  turbulent  wakes  might  not  be  homogeneous.  However, 
they  have  weaker  enstrophy.  Most  of  the  vorticity  variance  is  at  small  scales,  characteristic 
of  blue  spectra.  Therefore,  the  horizontal  homogeneity  assumption  should  be  satisfied 
approximately. 

We  computed  the  average  spanwise  and  vertical  components  of  enstrophy  in  each 
transect  when  the  ship  steamed  at  a  nearly  constant  velocity,  excluding  measurements 
when  the  ship  turned,  accelerated,  or  decelerated.  The  horizontal  and  vertical  components 
of  enstrophy  have  a  similar  magnitude  and  are  significantly  correlated  with  a  correlation 
coefficient  of  0.8  and  the  95%  significance  level  of  0.35  (Fig.  3a).  At  an  intermediate 
enstrophy  strenght,  the  vertical  component  is  slightly  greater  than  the  horizontal 
component.  The  estimates  of  total  enstrophy  agree  with  observed  e  (Fig.  3b),  as  suggested 
by  the  high  Reynolds  number  relation.  The  correlation  coefficient  between  the  logarithms 
of  £2  and  e  is  0.91,  with  the  95%  significance  level  of  0.35. 
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The  probability  density  functions  (pdf)  of  the  magnitudes  of  vorticity  components 
and  the  total  vorticity,  and  the  pdf  of  s1/2  have  an  exponential  form  [4],  i.e., 

p(X)  ~  exp(-pAy<j>. ),  (3) 

with  P  =  1.4,  where  X  represents  any  of  demeaned  \£z\ ,  \%h\,  \^\ ,  and  eYl ,  and  a\  represents 

the  standard  deviation  of  X.  The  pdfs  of  observed  £  and  vorticity  magnitudes  show  a  nearly 
identical  exponential  form  (Fig.  2c).  They  closely  follow  the  empirical  formula  (3).  These 
results  (Fig.  3)  confirm  the  quality  of  measurements  of  C,  and  e ,  and  form  the  basis  of  the 

following  analysis.  Furthermore,  the  close  agreement  concludes  that  the  empirical  pdf 
formula  derived  from  low  Reynolds  number  flows  exists  also  in  the  high  Reynolds  number 
field  flows. 

4.  Downstream  Variations  in  Turbulence  Kinetic  Energy  Dissipation  Rates 

A.  Observed  Turbulence  Kinetic  Energy  Dissipation  Rates  Along  the  Centerline 

We  compute  fin  1-s  intervals  and  examine  its  downstream  variation.  Because  of 
the  limited  spatial  resolution  of  our  measurements  to  resolve  the  cross-wake  and 
along-wake  structure  and  to  avoid  the  complexity  of  the  cross-wake  variation,  we  only 
investigate  the  downstream  variation  of  the  dissipation  rates  along  the  centerline  of  the 
wake.  The  Tacoma  Narrows  bends  at  -1800  m  downstream  from  the  pier  of  the  Tacoma 
Narrows  Bridge  (Fig.  la)  beyond  which  the  modulation  of  vorticity  and  turbulence  by  the 
lateral  boundary  cannot  be  ignored.  In  the  following  analysis,  we  discuss  only  observations 
within  1 800  m  downstream  of  the  pier. 
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Turbulence  kinetic  energy  dissipation  rates  s  observed  within  10  m  from  the 
centerline  of  the  wake  are  averaged,  typically  a  10-s  time  interval  with  the  ship’s  drift 
velocity  of  -4  kts  (Fig.  4).  Upstream  of  the  pier  the  background  turbulence  is  not  weak,  s 
~10“5  W  kg-1,  showing  the  turbulent  nature  of  the  flow  in  Tacoma  Narrows.  Such  an 
energetic  background  is  different  from  most  laboratory  experiments  in  which  a  weak 
turbulence  background  is  often  enforced. 

During  the  first  two  legs  of  the  experiment,  there  were  no  observations  at  the 
wake’s  centerline  within  400  m  downstream  of  the  pier.  At  -400  m  downstream  from  the 
pier,  e  is  20-50  times  of  that  observed  upstream  of  the  pier.  Further  downstream,  it  decays 
almost  monotonically  with  distance.  The  observed  e  has  a  similar  magnitude  during  the 
first  two  legs,  and  is  the  weakest  in  the  third  leg,  as  expected  in  the  weakening  ebb  tidal 
current. 


The  third  transect  shows  a  weak  value  of  s  at  300  m  downstream  and  the  strongest  e 
at  800  m  downstream.  Our  measurements  do  not  have  sufficient  spatial  and  temporal 
resolutions  to  investigate  this  departure  from  the  expected  monotonic  downstream  decay. 
The  turbulence  intermittency  and  unsteadiness  of  the  turbulent  wake  might  be  the  sources 
for  the  observed  discrepancy. 

The  dependence  on  the  tidal  forcing  is  better  illustrated  by  the  normalized  e  and  the 
normalized  downstream  distance,  scaled  respectively  as  e,  =  £/(ul,d-')  and  x.=x/d  (Fig. 
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4b).  The  choice  of  normalization  is  discussed  in  the  next  section.  This  normalization 
reduces  the  nearly  one-decade  scattering  of  e  to  within  a  factor  of  two  (Figs.  4a  and  4b). 

We  apply  two  types  of  curve  fittings,  power  law  and  exponential,  to  observations 
within  1 2  <x*  <60  to  yield  0.0  l(x.  /  xc)  3  and0.008e  (  *'_xc)/10  The  referenced  downstream 
scale  xc  =  12  is  chosen  where  the  observed  fis  the  strongest.  The  dimensionless  amplitude, 
0.01  and  0.008,  the  power  law  x“3 ,  and  the  exponential  decay  scale  10  d  are  parameters 
determined  by  the  least  squared  fit  to  observations.  The  obtained  -3  power  law  is  steeper 
than  that  predicted  by  the  self-preservation  similarity  in  the  turbulent  far  wake,  as  reviewed 
in  the  next  two  sections.  Sreenivasan  [5]  reports  that  the  decay  rate  of  a  temperature 
anomaly  in  a  heated  cylindrical  turbulent  wake  is  faster  in  the  intermediate  and  near  wake 
than  in  the  far  wake.  Our  observed  jc.“3  decay  rate  for  ein  the  intermediate  wake  does  not 
conflict  with  results  of  previous  studies. 

B.  Review  of  Self-preservation  Similarity  in  a  High  Reynolds  Number  Turbulent  Far 
Wake 

In  a  turbulent  wake,  mean  and  turbulence  properties  may  be  expressed  as  functions 
of  two  parameters  [3],  the  velocity  defect  at  the  centerline  of  the  wake,  us=ux-  u(y  =  0) , 

where  u(y= 0)  is  the  along-wake  velocity  at  the  centerline,  and  £,  the  half-wake  thickness 
defined  as  us  12  =  um  —  u(y  =  A) .  Both  €  and  us  vary  in  the  downstream  direction 

(x-component)  and  depend  on  the  obstacle  scale  d  and  the  ffee-stream  velocity  w® 
expressed  as  [6] 
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(4) 


X=Cd(^-) 


x-Xo  \l/2 


Us=CuU„(^, 


(5) 


where  x<>  is  the  virtual  origin.  Previous  studies  report  indefinite  values  of  the  virtual  origin, 
varying  from  -40  d  to  -125  d.  Some  argue  that  the  virtual  origin  depends  on  the  Reynolds 
number.  The  values  of  coefficients  Cx  and  Cu  are  estimated  as  0.2  and  1 .28  by  Browne  and 
Antonia  [6].  Similar  values  are  reported  elsewhere  [3]. 


The  along-wake  and  across-wake  self-preservation  similarity  of  velocity,  the 
dissipation  rates  of  turbulence  kinetic  energy,  and  enstrophy  are  described  as  [3] 

•  - =  /u(T.) 

us 

>  i  .  .  •  '  •  • 

*  =  ^/.00  =  C,  [u'jd\dAx-x0)Yay.) 

C2  =^ft(y-)  =  Ce  k  1  vd  \d  ~  x0  )]2  fs  (T.  )  » 


(6) 

(7) 

(8) 


where  y,  =  y"k  1  is  the  non-dimensional  cross-wake  coordinate,  and  Cc  =  C]  /  Cx. 

The  analytical  form  for  the  cross-wake  structure  fu(y.)  and  fc(y.)  have  been  discussed 
previously  [7,  8, 9].  This  is  beyond  the  scope  of  our  study  and  will  not  be  discussed  further. 


Previous  experiments  conclude  that  the  self-preservation  properties  establish  at  x  > 
1 00  d  in  the  far  wake  [10,  11].  Our  observations  are  taken  at  x  <  60  d,  in  the  near  and 
intermediate  wake  regions  and  the  observed  turbulent  wake  is  unsteady  due  to  the  time 
varying  tidal  current.  We  do  not  expect  self-preservation  similarity  to  apply  to  our 
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observations.  Nonetheless,  we  construct  a  pure  kinematic  model  on  the  basis  of  the 
self-preservation  similarity  and  compare  model  results  with  observations,  not  to  justify  the 
similarity  scaling,  but  to  emphasize  the  difference. 


C.  Model  Prediction  of  e 

Following  the  steady-state  self-preservation  similarity,  we  construct  a  similarity 
scaling  imposed  by  a  time  varying  ffee-stream  velocity  as 

em(x,y,t)  =  g(u,(x,t)Nx,t))fs(y.),  (8) 

where  em  is  the  model  dissipation  rate,  g  is  a  function  of  time  varying  velocity 
scale us{x,t) and  the  wake  scale X(x,t) .  We  assume  that  us(x,t )  and  A(x,0  scale  by  the 

time  varying  free-stream  velocity  u"  (t),  and  the  obstacle  scale  d.  The  proposed  scaling  for 
e  is  described  as 


1-2 


i 

f  u”(t‘)dt'—x0 


/.O'.)  =  #,(*,<)/,  O’.). 


(9) 


The  downstream  distance  from  the  pier  is  x  =  fir”(r,)rfr’.  where  t0  is  the  time  when  the 

*0 

turbulence  is  generated  at  the  pier.  Turbulence  is  advected  downstream  by  the  decaying 
tidal  current. 


Turbulence  kinetic  energy  dissipation  rate  predicted  by  the  above  equation  for 
*o  =  — 40c/  is  shown  in  Fig.  5.  The  model  results  of  e  at  the  centerline  are  compared  with 
our  observations  (Fig.  4).  Model  results  decrease  by  nearly  one  decade  during  the  3-h 
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experiment  period,  similar  to  the  range  of  observed  variation  (Fig.  4a),  and  have  a 
magnitude  comparable  to  our  observations.  However,  the  downstream  variation  of  em  is 

much  slower  than  the  observed  x'3  variation.  Model  em  varies  as  x~2  in  the  far  field,  x  »  xo, 
far  beyond  our  observation  regime.  Model  results  using  different  virtual  origins  reveal  the 
similar  slower  downstream  decay  rate  than  the  observed  x'3  decay  rate. 

The  observed  x'3  cannot  be  explained  by  the  time  varying  free-stream  velocity.  We 
propose  that  x“3  is  the  true  decay  rate  of  dissipation  rates  of  turbulence  kinetic  energy  in 
the  intermediate  wake.  The  faster  decay  rate  of  turbulence  is  further  confirmed  by  the 
variation  of  the  observed  turbulence  enstrophy. 

5.  Enstrophy  Properties 

A.  Observations  of  Enstrophy  Along  Centerline 

The  downstream  variations  of  the  vertical  and  horizontal  components  of  enstrophy 
along  the  centerline  of  the  turbulent  wake,  within  10  m,  are  shown  in  Figs.  6  and  7.  The 

_2 

background  value  of  both  components  of  enstrophy,  upstream  from  the  pier,  is  about  1  s  . 
At  -300  m  downstream  from  the  pier,  both  components  of  enstrophy  increase  by  more  than 
two  decades  above  the  upstream  intensity.  They  decrease  downstream  nearly 
monotonically  and  by  more  than  one  decade  at  1600  m  from  the  pier.  Both  components  of 
enstrophy  during  leg  2  are  stronger  than  those  during  leg  3,  as  expected  by  the  tidal  forcing. 
Unexpectedly,  both  components  of  enstrophy  during  leg  1  are  weaker  than  those  during  leg 
2.  We  do  not  have  an  explanation  but  attribute  this  to  turbulence  randomness. 
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We  normalize  the  observed  enstrophy  by  u\ldv,  in  accordance  to  the  scaling  for 
e.  This  normalization  does  not  collapse  the  scattering  significantly,  partially  due  to  the 
unexpected  low  enstrophy  during  leg  1.  The  power  law  fit  yields  oc  X;2A  and  oc  X;21  for 
the  vertical  and  the  horizontal  components  of  enstrophy,  respectively,  in  close  agreement 
with  the  downstream  decay  rate  of  e.  The  exponential  fit  yields  the  decaying  scale  of ~10  d, 
again  consistent  with  that  of  e. 

B.  Turbulence  Intermittency 

The  spatial  spottiness  of  the  fine  scale  structure,  the  so-called  internal  intermittency, 
is  a  well-known  feature  of  turbulence.  Previous  studies  often  quantify  the  intermittency 
using  the  decay  rate  of  the  autocorrelation  of  turbulence  properties  with  increasing  lags. 
Shafi  et  al.  [12]  investigate  the  intermittency  within  a  turbulent  wake  using  the  structure  of 
the  autocorrelation  of  enstrophy,  i.e., 

^(/?)  =  <4'2(rK2(r  +  /?)>/^J  ~(L/Ry,  (10) 

where  the  exponent//  is  the  intermittency,  R  is  the  spatial  lag,  and  L  is  the  turbulence 
integral  length  scale.  Previous  studies  report  the  intermittency  //  to  be  about  0.2  in  low 
Reynolds  number  flows  [12,  13]. 

We  compute  the  autocorrelation  coefficients  using  observed  enstrophy  within  10  m 
from  the  centerline  and  within  10  and  60  d  downstream  from  the  bridge  pier.  Observations 
taken  at  irregular  spatial  intervals  are  first  interpolated  to  a  constant  spatial  interval  ~  0. 1  m 
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before  computing  the  autocorrelation  coefficients  y  2  .  Values  ofy(2 ,  averaged  within 

prescribed  intervals  of  spatial  lags,  of  both  enstrophy  components  exhibit  a  monotonic 
decaying  power  exponent  close  to  those  reported  previously  [12].  During  our 
experiment,  e  varies  between  10"*  and  10"3  W  kg"1,  and  77  varies  from  0.03  cm  to  0.1  cm. 
The  autocorrelation  coefficient  decreases  at  a  nearly  constant  rate  for  lags  between  0.3  m 
and  4  m,  corresponding  to  300  <  RJr\  <  1 04.  We  interpret  this  as  being  part  of  the  turbulence 
inertial  subrange.  Outside  of  this  range,  the  autocorrelation  coefficient  exhibits  different 
slopes  with  respect  to  the  spatial  lags.  The  power  exponent  decay  rates,  i.e.,  the 
intermittency  p.  for  v  ,  and  v  ,  within  0.3-4-m  scales,  are  0.2  and  0.22,  respectively, 

computed  by  fitting  logarithmic  values  of  y  t  and  R  using  the  least  squared  method.  These 

results  are  in  near  agreement  to  those  reported  in  previous  studies  of  turbulent  wakes  and 
jets  [12]. 

6.  Summary 

Dissipation  rates  of  turbulence  kinetic  energy  and  the  vertical  and  horizontal 
components  of  enstrophy  were  measured  in  a  turbulent  wake  behind  a  bridge  pier.  The 
turbulent  wake  is  generated  by  the  tidal  current.  Our  observations  confirm  the  relation 
between  the  enstrophy  and  dissipation  rates  of  turbulence  kinetic  energy  in  high  Reynolds 
number  flows,  and  confirm  the  empirical  probability  density  functions  reported  in  previous 
laboratory  studies. 
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Dissipation  rates  of  turbulence  kinetic  energy  and  two  components  of  enstrophy 
near  the  centerline  exhibit  a  similar  downstream  decay  rate,  one  faster  than  that  predicted 
by  the  self-preservation  similarity  in  the  far  wake.  We  propose  that  the  larger  decay  rate 
exists  in  the  intermediate  wake. 

The  turbulence  intermittency  is  estimated  using  the  autocorrelation  function  of  the 
two  components  of  enstrophy.  Autocorrelation  coefficients  of  enstrophy  exhibit  distinct 
slopes  in  different  ranges  of  spatial  lags.  In  the  scale  much  greater  than  the  Kolmogorov 
length  scale  and  smaller  than  4  m,  identified  as  a  part  of  the  inertial  subrange,  the 
intermittency  is  close  to  0.2,  consistent  with  laboratory  experimental  results. 
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FIG.  1 .  (a)  Bathymetry  of  Tacoma  Narrows  and  ship  tracks.  Land  is  shades  of  gray;  water 
bathymetry  is  all  other  colors  colors.  Thin  contours  are  1 0-m  isobaths;  the  thick  black 
contour  is  50  m  depth.  The  red,  yellow,  and  magenta  zigzag  curves  are  ship  tracks  during 
the  first,  second,  and  third  surveys  (as  labeled),  respectively,  downstream  of  a  pier  of  the 
Tacoma  Narrows  Bridge  (shown  as  a  black  line  with  white  edges).  Values  inside 
parentheses  indicate  the  total  sampling  period  in  seconds.  The  thick  yellow  arrow 
illustrates  the  direction  of  the  ebb  tide.  The  two  perpendicular  red  arrows  indicate  the 
coordinate  system  in  the  downstream  and  cross  stream  directions,  (b)  Photograph  of 
vortices  behind  the  pier  of  the  Tacoma  Narrows  Bridge  taken  soon  after  the  bridge 
collapsed  in  1940.  (c)  Sketch  of  Electro-Magnetic  Vorticity  Meter  (EMVM).  Two  vorticity 
meters  (VMs)  are  mounted  on  the  side  and  the  bottom,  respectively.  A  microstructure  shear 
probe  is  mounted  in  the  front  of  the  EMVM.  CTD  sensors  and  a  pressure  sensor  are 
mounted  on  the  side  of  the  EMVM.  Other  measurements  include  compass,  pitch,  pitch  rate, 
and  yaw  rate  of  the  EMVM  to  calibrate  the  instrument  motion. 
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FIG.  2.  Temporal  and  transverse  variations  of  the  free-stream  velocity.  Filled  circles  are 
averages  of  observed  streamwise  velocity  outside  of  the  turbulent  wake.  Grey  shadings 
represent  the  transverse  distance  from  the  centerline  of  the  wake.  The  solid  curve  is  the  fit 
of  observations  to  a  model  current  field  consisting  of  a  transverse  shear,  a  mean  current, 
and  a  semidiurnal  tidal  current.  Dashed  curves  represent  the  95%  confidence  interval  of  the 
model  fit. 
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FIG.  3.  (a)  Comparison  of  the  horizontal  and  vertical  components  of  vorticity  variances 
measured  from  the  EMVM,  (b)  comparison  of  estimates  of  total  enstrophy  and  turbulence 
kinetic  energy  dissipation  rates,  and  (c)  probability  density  functions  of  the  normalized 
squared  root  of  e  (solid  dots),  demeaned,  normalized  magnitudes  of  horizontal  vorticity 
(open  circle),  of  vertical  vorticity  £z(open  square),  and  of  total  vorticity  C,  (black  diamond). 
The  black  curve  is  the  empirical  curve  [4],  The  variable  a  represents  demeaned 
£l/2,  |^h|,  \Q,  and  ICI>  and  aa  represents  its  standard  deviation. 
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FIG.  4.  (a)  Downstream  variation  of  turbulence  kinetic  energy  dissipation  rates  along  the 
centerline  of  the  turbulent  wake  (solid  lines  with  symbols),  and  (b)  the  non-dimensional 
form  of  (a)  where  the  downstream  distance  is  normalized  by  the  scale  of  the  pier  d,  and  the 
turbulence  kinetic  energy  dissipation  rates  normalized  by  w*  Id  where  ux  is  the 
free-stream  velocity.  The  black,  dark  grey,  and  light  grey  curves  in  (a)  are  model  results 
for  the  three  transects  obtained  by  assuming  the  property  of  self-preservation  similarity. 
The  black  and  grey  solid  curves  in  (b)  are  power  law  fit  and  exponential  law  fit  to 
observations. 
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FIG.  5.  Model  spatial  and  temporal  variations  of  s  induced  by  the  ebbing  tidal  current 
interacting  with  the  bridge  pier.  The  tidal  current  model  shown  in  Fig.  2  is  used.  Three 
thick  dashed  curves  represent  survey  measurement  tracks. 
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FIG.  6.  (a)  Downstream  variation  of  the  vertical  component  of  enstrophy  along  the 
centerline  of  the  turbulent  wake  (solid  lines  with  symbols),  and  (b)  the  non-dimensional 
form  of  (a)  where  the  downstream  distance  is  normalized  by  the  scale  of  the  pier  d  and  the 
enstrophy  is  normalized  byw*  Ido.  The  black  and  grey  solid  curves  in  (b)  are  the  power 
law  fit  and  the  exponential  law  fit,  respectively,  to  observations  in  10  <x/d<  60. 
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FIG.  7.  (a)  Downstream  variation  of  the  horizontal  component  of  enstrophy  along  the 
centerline  of  the  turbulent  wake  (solid  lines  with  symbols),  and  (b)  the  non-dimensional 
form  of  (a)  where  the  downstream  distance  is  normalized  by  the  scale  of  the  pier  d,  and  the 
enstrophy  is  normalized  by  u^/du.  The  black  and  grey  solid  curves  in  (b)  are  the  power 
law  fit  and  the  exponential  law  fit,  respectively,  to  observations  in  10  <x/d  <  60. 
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FIG.  8.  Autocorrelation  coefficients  of  vertical  enstrophy  £]  (solid  curve),  and  of 
horizontal  enstrophy  (dashed  curve).  Shadings  are  their  95%  confidence  intervals.  The 
regression  power-law  fits  in  0.3  <  R  <  4  m,  where  R  is  the  lag,  are  labeled. 
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Table  1.  Measurements  and  sensors  mounted  on  the  Electro-Magnetic  Vorticity  Meter. 


Measurement 

Manufacturer 

Sample  rate  (Hz) 

Resolution 

Vorticity  (two) 

APL-UW 

20 

3  xlO-4  s l 

Velocity  (two) 

APL-UW 

20 

10-5  m  s-1 

TKE  dissipation 

APL-UW 

400 

10-8  W  kg-1 

Temperature 

Sea-Bird 

20 

0.5  *10~3K 

Conductivity 

Sea-Bird 

20 

0.5  *10”4  S  m”1 

Pressure 

Paroscientific 

20 

0.1  db 

Magnetic  heading 

KVH,  Inc. 

1 

0.1° 

Pitch 

IC  sensors 

20 

0.035° 

Roll 

IC  sensors 

20 

0.035° 

Pitch  rate 

Humphrey 

20 

10"3°  s"1 

Yaw  rate 

Humphrey 

20 

10"3°  s”1 

Vessel  Position 

Northstar 

0.5 

3-5  m 
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